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ABSTRACT 

The soft X-ray background in the 0.1- 1.0 keV band is known to be produced by at least 
three sources; the Local Hot Bubble (LHB), the extragalactic power law (EPL), and a seemingly 
galactic component that lies outside the bulk of the absorption that is due to the ISM of the 
galactic disk. This last component, which we call the Trans- Absorption Emission (TAE), has 
been modeled by a number of groups who have derived disparate measures of its temperature. 
The differences have arisen from differing assumptions about the structure of the emitting gas and 
unrecognized methodological difficulties. In particular, spectral fitting methods do not uniquely 
separate the TAE from the foreground emission that is due the LHB. This “degeneracy” can be 
resolved using the angular variation of the absorption of the TAE. 

We show that the TAE cannot be characterized by a single thermal component; no single- 
component model can be consistent with both the spectral energy distribution of the TAE emis- 
sion and the angular variation due to absorption by the galactic disk. We use the angular 
anticorrelation of the ROSAT All-Sky Survey with the galactic absorption to separate local from 
distant emission components, and to fit the spectral energy distribution of the resulting distant 
emission. We find that the emission is best described by a two-thermal-component model with 
log 7s = 6.06 ^q \ j and logT# = 6A2to'\2- This two-thermal-component TAE fits the ROSAT 
spectral energy distribution significantly better than single-component models, and is consistent 
with both angular variation and spectral constraints. 

Subject headings : galaxies: Milky Way — interstellar: matter — X-rays: galaxies, general 


1. Introduction 

The canonical model of the soft X-ray back- 
ground (0.1- 1.0 keV) posits a Local Cavity (in the 
H I of the galactic disk) that is filled, at least in 
part, by hot gas emitting mostly in the \ keV 
band. This region of X-ray emitting gas is known 
as the Local Hot Bubble (LHB). The LHB emis- 
sion is absorbed by the nearby very tenuous Local 
Interstellar Cloud (LIC), otherwise known as the 
Local Fluff. Beyond the Local Cavity there are 
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thicker layers of absorbing gas which partially ab- 
sorb emission from an isotropic extragalactic com- 
ponent and emission from at least one other com- 
ponent. This latter component we will call the 
“trans-absorption emission” (TAE) or the “dis- 
tant component”. The local component and the 
distant component are usually modeled as ther- 
mal plasmas. The isotropic extragalactic compo- 
nent, which is most clearly seen at E > 1 keV, 
and is presumably due to the accumulation of un- 
resolved point sources, is modeled as a power-law, 
and will be referred to as the Extragalactic Power 
Law (EPL). The distant component is required to 
contribute to the observed excess over the extrap- 
olation of the EPL in both the \ keV and 7 keV 
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bands. 

Table 1 contains a compendium of measures of 
the temperature of the thermal components us- 
ing various techniques. The ROSAT measures 
alone show a wide dispersion. In particular, PSPC 
spectral fitting produces lower LHB temperatures 
and higher halo temperatures than do early broad- 
band measures and PSPC broad-band deabsorp- 
tion measures. In order to understand the source 
of this difference, the next section assumes the 
canonical model and explores the limitations of 
the methods used for determining the parameters 
of that model. 

We find that although a portion of the differ- 
ence is due to inadequacies in the spectral fitting 
method, the bulk of the difference is due to as- 
sumptions about the nature of the distant compo- 
nent. Most investigators have assumed that the 
TAE is due to a single-temperature galactic halo. 
We find that none of the existing single-component 
models satisfies all of the constraints placed upon 
the data by the angular variation of the absorption 
of the TAE. 

We attempt to determine the thermal structure 
of the TAE by combining the all-sky deabsorption 
method of Snowden et al. (1998) with six-band 
Spectral Energy Distribution (SED) fitting. Sec- 
tion 3 describes the assumptions, input data, and 
limitations of our method. We find that although 
the data require more than one thermal compo- 
nent to describe the TAE, the data are sufficient 
to characterize only two components, though a 
greater number may exist, and the temperature of 
these components may vary with position. Section 
4 describes the results and finds the temperatures 
of the two thermal components. Section 5 summa- 
rizes our arguments and discusses the implications 
of a two-component TAE. 

2. Analysis of Previous Results 

Using the canonical model, several efforts have 
been made to determine the spectrum of the soft 
X-ray background. We divide these methods into 
two groups; spectral fitting of individual PSPC 
pointings, and all-sky deabsorption of broad-band 
images. In both of these cases, the data were fit 
to the function 

Io = I L e- {Nl)<r{v ’ N >' S] 


+/£)e - *' Vfc) ‘ T * ,/ '' v ‘' ,s ^ (1) 

A r b )*(*/, .^,5) 

where Io is the observed X-ray flux. I l is the lo- 
cal (LHB) X-ray emission which may be absorbed 
by an intervening layer of column A’/, Io is the 
distant emission (TAE) and Ie is the contribu- 
tion due to the EPL. Both the distant and extra- 
galactic components are absorbed by a column AV 
The quantity a is the absorption cross-section. In 
the case of spectral fitting, cr is dependent only 
upon energy (//), but for broad-band deabsorption 
we use the effective cross-section for a given band 
which is dependent both upon the magnitude of 
the absorbing column (A r ), and upon the spectral 
shape of the background emission (d>). 

2.1. Spectral Fitting 

Several groups (e.g., Chen et al. 1997 and 
Miyaji et al. 1998) have attempted to fit the spec- 
trum of the X-ray background as observed in in- 
dividual PSPC pointings. Attention has focused 
upon determining the index of the EPL; determi- 
nation of the temperatures and normalizations of 
the two thermal components has been incidental. 
Due to this focus, the ROSAT PSPC data have 
usually been supplemented wdth A SC A observa- 
tions of the same field in order to gain a greater 
lever arm for fitting the EPL. Fitting the thermal 
components is more problematic. 

The thermal components occupy only a very 
limited spectral region in the ROSAT range; the 
emission lies mostly within 0. 1-0.4 ke\ r , which is 
completely outside of the A SC A range. At these 
energies, absorption, and the uncertainty due to 
its correction, is much greater. Also greater are 
the non-cosmic foregrounds; the Long Term En- 
hancement (LTE) and the Scattered Solar X-rays 
(SSX). The LTE is a particular problem for PSPC 
pointed observations as adequate removal requires 
a light curve extending over (at minimum) several 
weeks to months (Snow’den et al. 1995, 1994b). 
Much of the disagreement between different spec- 
tral fits has been attributed to inadequate non- 
cosmic background removal (Miyaji et al. 1998). 

Inadequate spectral resolution is also a prob- 
lem. There are only about five independent spec- 
tral resolution elements in the ROSAT PSPC en- 
ergy band, and only three spanning the energies 
emitted by the thermal components. Thus, even 


2 


II 


Table 1 
Compendium 


Source 

t l 


T d i 


Td 2 


keV 

logT 

keV 

logj 

keV log J 


Pre-ROSAT Results 

Block et al. (1986) 

0.077-0.097 

5.95-6.05 




Juda et a!.(1991) 

0.054-0.136 

5. 8-6. 2 




Snowden et al. (1990) 

0.086 

6.0 




ROSAT Era Results 

Chen, Fabian k Gendreau(1997) 

0.05 

5.76 

0.12 

6.14 


Kerp(1994) 

n i +0.01 6 

0.061. _q 015 

5.85 ± 0.1 

H 1 79+0.045 
U ' A ' z -0.035 

6.3 ±0.1 


Miyaji et al.(1998) 

.057-. 09 

5.82-6.02 

0.142-0.148 

6.217-6.233 


Peitz et al.(1998) b 



0.135 

6.195±0.05 b 


Sidher(1996) c 



0.137(0.216) 

6.2 (6.4) c 


Snowden et al.(1998) d 

0 1 1 Q+° 034 

U. I iy_o.012 

6.14t°;J^ 

0 106 +o 022 
U.lUO_o .018 

nq+ 0 08 e 
D,uy -0 .08 


Snowden et al.(2000) h 

0 1 1 Q +0 034 
uliy -0.024 

6.Ut° 0 \ l 0 efg 

0 100+° 033 
U.iUU_o 030 

6 07 +0 12e 
o.u / _ 0 .16 


ROSAT Results - This Work 

Hyades 

0. 103^0 QQ2 

6.08±° gp 




Galactic Plane 

0.108^007 

pt 1 1 +0.02ei 
O.i 1_ 0 03 




Section 3.4 

O.lOQtHta 

e.ioi^p/ 

0.098°« 3 o 5 31 

6.06ppp 

O 997+0.086 yj9+0.14e 

l — 0.055 °-^ z -0.12 

Section 4.1.4 

0.108 ±0.04 

6.10 ±0.02 

0.087 ± 0.003 

6.00 ±0.02 

0.180 ±0.008 6.32 ±0.02 

Section 4.1.4 

0.127 ±0.010 

6.17 ±0.03 

0.313 ± 0.055 

6.56 ± 0.07 


Post-ROSAT Results 

Parmar et al.(1999) k 

0 1Q7+ 0011 
U.io / _o .022 

a O+0.03 
°' Z -0.07 

0 7 +0 4 

U * ' -0.3 

6-91^ k 



a BB * Broad-Band Ratios, SF - Spectral Fitting, DA - De-Absorption 

b Temperature derived from the Mewe-Kaastra (Mewe et al. 1985) code. A fit of a Raymond k Smith (1977) spectra to 
Kaastra spectrum absorbed by 1.5 x 10 20 cm -2 and folded through the ROSAT PSPC response yielded a Raymond & Sm 
temperature of logT = 6.18. If a Mewe-Kaastra- Leidahl spectrum were used rather than a Mewe-Kaastra spectrum, the R 
Smith plasma temperature would be logT = 6.17. 

c Temperatures derived from the Landini &: Monsignori-Fossi (Landini k Monsignori-Fossi 1990) code. Sidher et al. note tha 
k Smith code would produce a log To == 6.4. 

d A second, hotter component is implied, but not explicitly characterized in this paper. The values given here are not quit* 
given in the original paper. Those values had been calculated from band ratios using an older detector response matrix; wc 
the current detector response matrix. 

e Cited range is the FWHM of the distribution. These distributions are dominated by fit and/or data uncertainties; the c 
of temperatures is very narrow. See text §4 for a fuller explanation. 

f The upper limit is not given by the upper FWHM of the distribution but by the highest R2/R1 ratio permitted by an i 
Raymond k Smith spectra. 

g These LHB temperatures have not taken into account absorption within the LHB. Assuming an absorbing column ~ 4 x 
the temperatures drop by AlogT ~ 0.01. 3 

h This paper assumed a second, hotter component with \ogTo 2 ~ 6.4. 

'This temperature was derived in the anti-center. Snowden et al. (1990b) showed that the temperature of the LHB dej 
galactic longitude; with galactic center directions having slightly higher temperatures than the galactic anti-center. 

k Result obtained with the BeppoSAX LECS. The temperatures were derived using the Mewe-Kaastra-Liehdal code. T1 
should not be compared to the current work because the Parmar et al. (1999) data selection does not well sample the gi 
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Fig. 1. — A comparison of the spectra of ther- 
mal components with the available independent 
spectral elements for the ROSA TPSPC. Top: the 
ROSAT PSPC band response functions. Below : 
Three thermal spectra spanning the range of tem- 
peratures found in the literature. The vertical 
lines show' the divisions between independent spec- 
tral resolution elements, i.e., elements with non- 
overlapping FWHM. The bulk of the emission 
from galactic thermal components falls almost en- 
tirely within two independent spectral resolution 
elements. 

given some a priori model of the EPL, the four 
fit parameters necessary to describe two thermal 
components are under-determined. In practice, 
the situation is even worse. Figure 1 shows the 
Raymond & Smith (1977) spectra for two repre- 
sentative thermal components, and a maximal set 
of independent resolution elements. (As with all 
future references to the Raymond &: Smith spec- 
tra in this paper, except in the context of XSPEC, 
we are using the 1991 release with cosmic abun- 
dances.) The bulk of the emission is contained in 
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Fig. 2. — The degeneracy in the thermal com- 
ponent fit parameters. We have used a “mean” 
high galactic latitude background spectrum wdth a 
mean absorbing column (Nt,) of 1.21 x 10 20 cm" 2 
(and Nj =0), have assumed a Chen, Fabian & 
Gendreau (1997) EPL (their Model A), and have 
assumed a single-temperature distant component. 
Top: The range of normalizations required for the 
local emission component, in units of 1000 n e n p L 
( n e and Tih in cm -3 and L in pc), as a function 
of the local and distant component temperatures. 
Middle: The same as above for the distant com- 
ponent. Bottom: The region of parameter space 
allowed by the R4, R5, and R6 band fluxes, allow'- 
ing a ±5x 10“ 6 counts s~ J arcmin“ 2 uncertainty 
in the fluxes. Solid: R4, Dotted: R5, and Dashed: 
R6. Results are similar for other EPL models, 
though no EPL model can force the three regions 
to have significant overlap. 

— 2 resolution elements. 

At the risk of appearing to labour the point, 
consider the following. The six standard ROSAT 
PSPC bands R1-R7 (the R1 and R2 bands form- 
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ing the | keV band, the R4 and R5 bands form- 
ing the | keV band, and excluding R3 as usual, 
Snowden et al. 1995) contain all of the information 
in the PSPC spectrum. The bulk of the informa- 
tion about the thermal components is contained in 
the R1 and R2 bands. Using only the R1 and R2 
bands, given the absorbing column and the EPL, 
one can choose (within reason) any pair of local 
and distant temperatures, Ti and To, and find 
some pair of normalizations Nl and No that pro- 
vide the observed R1 and R2 band values. Figure 2 
shows the ranges of normalizations required over 
a reasonable range of Ti and To 4 > 


Also shown in Figure 2 is the region allowed by 
the R4, R5, and R6 band values. The horizontal 
branch of the allowed parameter space represents 
a “minimal LHB” model, where the £ keV flux is 
dominated by the distant emission; the LHB com- 
ponent is then poorly constrained and can assume 
a broad range of values. The vertical branch of 
the allowed parameter space represents a “maxi- 
mal LHB” model, where the ^ keV flux is dom- 
inated by the LHB; the distant emission is then 
poorly constrained and can assume a broad range 
of values. It should also be noted that the vertical 
branch falls at the Ti at which the R2/R1 ratio for 
the local component is a maximum (see Figure 3.) 

Most of the published parameters of the local 
and distant thermal components fall in the region 
allowed by the R4 flux. 


The “degeneracy” in the fit parameters can be 
broken with the addition of angular information; 
if Nb varies over scales smaller than the angular 
scales of the variation in the emission, then the lo- 
cal and distant normalizations can be determined, 
and from them, the temperatures. It is not our 
desire to single out any particular practitioner of 
spectral fitting for criticism, especially since the 
thermal parameters were not their prime objec- 
tive. However, one can easily show that the ther- 
mal components fit, for example , by Chen et al. 
(1997) are incorrect. Taking a small strip of the 
RASS in the southern galactic hemisphere that is 
free of strong ~ keV structures (t = 165°, —55° > 
b > —75°), and assuming that the temperatures 
of the local and distant emission do not vary with 


4 The data used in the construction of Figure 2 were regions 
with 6 > 45°, 45° < t < 270°, and 1.1 x 10 20 < N H < 
1.3 x 10 20 cm” 2 . 
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Fig. 3. — The R2/R1 ratio as a function of tem- 
perature. The top curve is the R2/R1 ratio for the 
distant component after absorption by a “mean” 
high galactic latitude column of 1.21 x 10 20 cm’ 2 ; 
the bottom curve is the ratio for the local com- 
ponent. The horizontal line is the ratio for the 
“mean” high galactic latitude background. Note 
that since the local component has an R2/R1 ra- 
tio less than the observed R2/R1 ratio in all cases, 
the distant component must have a temperature 
greater than logT* ~ 5.95. 


position, one can solve for the normalizations as 
a function of position. A plot of the normaliza- 
tion as a function of absorbing column density, 
Figure 4, shows a strong correlation between the 
absorbing column density and the normalization 
of the distant component. One should not expect 
the emission from beyond the absorbing layer to be 
correlated with the intervening column density of 
H I. Therefore, the normalization must be wrong. 

Indeed, repeating the exercise using the param- 
eters determined by Snowden et al. (1998), shows 
no correlation between the absorbing column den- 
sity and the normalization for the distant com- 
ponent. The uncertainties for the Snowden et al. 
(1998) parameters are larger than for the Chen 
et al. (1997) parameters because a softer distant 
component, absorbed by the intervening column, 
will require a greater change in No before effecting 
the observed flux. 
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Fig. 4. — Checking spectral fit parameters using 
absorption. Top: We have used the temperatures 
of Chen, Fabian & Gendreau (1997) (log 77, = 5.76 
and log To = 6.14) to demonstrate how the an- 
gular variation of the absorption can be used to 
check the temperature and normalizations derived 
from spectral fitting. Since one does not expect 
the true magnitude of the TAE to be correlated 
with the absorbing column, the fact that the de- 
rived TAE normalization is correlated with Nh in- 
dicates that the assumed component temperatures 
are incorrect. Bottom: We have used the temper- 
atures of Snowden et al. (1998) (log 7/, =6.12 and 
log To — 6.09). There is no visible correlation be- 
tween the TAE normalization and the absorbing 
column. 


Consider a spatially invariant background ab- 
sorbed by a spatially varying column of H I, to 
which is added a spatially invariant foreground. 
If one assumes that the background component 
has a higher temperature than the true back- 
ground, then the assumed effective absorption 
cross-sections will be lower than the true effective 
absorption cross-sections. Thus the expected vari- 
ation due to absorption (in the value of the trans- 
mitted flux) will be smaller than the true varia- 
tion, and that part of the true variation that is 
not accounted for by the assumed variation will 
be attributed (by whatever fitting mechanism one 
is using) to a variation in the strength of the back- 
ground component. Thus, if the assumed temper- 
ature of the distant component is wrong, or if the 
magnitude of the distant component is wrong for 
any reason, the fitted normalization will be corre- 
lated with the absorbing column. 

2.2. All-Sky Deabsorption 

The second approach to determining the X- 
ray background spectrum utilizes the galactic ab- 
sorption as a tool to separate the local emission 
from the distant emission. One assumes that at 
some angular scale the foreground and background 
emissions are constant at some II and Id, and 
that the variation in observed X-ray flux is due 
to absorption by the intervening column. Given a 
tracer of the absorbing column (e.g., either an H I 
map or an IRAS 100 pm map) one fits the relation 
in Equation 1 to find Ii and Id ■ 

Since we have very little knowledge of the distri- 
bution of N f , the absorption due to clouds within 
the LHB, one fits V L and Id in 

Io - I E e~ f - Nk ^ a ^‘'' Nb ' S ^ = I' L + 

(2) 

where 

r L zzI L e- (N ' )a(v ' N >' S) . (3) 

When using the I f L determined from several bands 
to solve for the temperature and emission measure 
of the local component, a uniform Nj is generally 
assumed (see § 4.2). Nj is small so, at worst, it 
can contribute a few percent error. 

2.2.1. \ keV (The R12 Band) 

Snowden et al. (1998) used the single band 
RASS maps of Snowden et al. (1997) to separate 
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the \ keV (R12 band) emission into foreground 
and background components over the entire sky. 
Although these maps have a mean exposure much 
lower than pointed observations, they have the ad- 
vantage of having the best possible LTE and SSX 
removal (Snowden et al. 1995), far better than 
can be accomplished for any single typical PSPC 
pointing. 

Snowden et al. (1998) assumed that both the 
local emission and the distant emission were flat 
on ~ 4?8 scales, and they fit Equation 1 assuming 
Nf = 0. The fit for each 4°8 x 4?8 region was then 
applied to center ninth of that region. They used 
the IRAS 100 /xm map of Schlegel et al. (1997) as 
a tracer of the absorbing column. This map was 
scaled by the mean value of H l/I\oo where the 21 
cm data (-74 km s“* < V < 24 km s -1 ) was taken 
from the Dwingeloo survey of Hartmann & Burton 
(1997). The EPL was assumed to have a T of 1.96 
(as per Hasinger et al. 1993) but was normalized to 
be consistent with the observed upper limit to the 
extragalactic background in the R12 band (Barber 
& Warwick 1994; Snowden & Pietsch 1995; Barber 
et al. 1996; Cui et al. 1996); although an inconsis- 
tent procedure, the inconsistency is inconsequen- 
tial, as shown by Snowden et al. (2000). The de- 
composition of the emission into foreground and 
background components was done separately for 
the R1 and R2 bands, and the R2/R1 band ratio 
was then used to determine the temperatures of 
the local and distant components. 

This method works well if the emission and ab- 
sorption is well described by the canonical model. 
If, however, the emission and the absorption are 
interleaved, then significant overestimates of the 
distant component will occur. Conversely, if the 
absorption tracer underestimates the total absorb- 
ing column, then the distant emission will be un- 
derestimated. Both problems pertain to regions 
with high and intermediate velocity complexes. 
(See Herbstmeier et al. 1995 and Kerp et al. 1999 
for preliminary studies of this problem.) 

As can be seen in Snowden et al. (1998), Figure 
8, (and again in our Figure 9) the distant emis- 
sion appears to be strongly mottled, even away 
from the problematic high velocity clouds, and 
there is evidence of temperature variation. Sur- 
prisingly, none of the structure observed in the 
t keV TAE is seen in the f keV band, which is 
smooth (see Figure 4 in Snowden et al. (1997) or 


our Figure 9). This difference led Snowden et al. 
(1998) to conclude that the component responsible 
for the structure in the \ keV band must not have 
significant emission in the | keV band, and must 
be different from the component producing the 
very smooth | keV emission. Accordingly, they 
felt justified in measuring the temperature of the 
component producing the | keV structure from 
the R2/R1 band ratio alone, leaving the measure 
of the hotter component to groups concentrating 
on the | keV band. 

Note that almost all other measures of the tem- 
perature of the TAE assumed a single-temperature 
component, and therefore measure a significantly 
higher temperature than Snowden et al. (1998). In 
this case, visible structure in the \ keV (R12) band 
but relative smoothness in the | keV (R45) band 
is understandable; for logT = 6.14, R45/R12 
~ 0.026, so much of the structure visible in the 
7 keV band will be reduced to the level of the 
noise in the 7 keV band. 

Since the effective cross-section for a given en- 
ergy band depends upon both the column density 
of the absorbing gas and the spectral shape of the 
background emission, one must assume a spectral 
shape (or temperature) for the distant emission 
before separating the observed emission into local 
and distant components. After the local/distant 
separation one can take the band ratios of the dis- 
tant component to determine a new spectral shape 
for the distant emission, and repeat the process. 

Did the assumption of the temperature of the 
distant emission bias the subsequent measure of 
the temperature of the distant emission? For a 
range of local and distant temperatures, we as- 
sumed a fixed single-temperature distant compo- 
nent, and determined the normalizations required 
to produce the total observed R1 and R2 band 
emission. Using those temperatures and normal- 
izations we then created vectors of flux vs. ab- 
sorbing column for both bands. Using a method 
similar to that used in Snowden et al. (1998), 
those vectors were then fit to derive the local 
and distant emission values, assuming the cross- 
sections appropriate for a distant component with 
log Id = 6.09. We used the R2/R1 band ratio 
from the resulting components to compute a “de- 
rived” temperature. Figure 5 shows the “derived” 
temperature for a range of true (or input) local 
and distant temperatures. One can see that for 



many temperatures, the Snowden et al. (1998) 
method (without iteration) could produce a bi- 
ased measure of the temperatures. However, for 
temperatures near the assumed temperature, the 
bias is small and does not preclude finding the 
true temperature through iteration. For all rea- 
sonable cases, the biases are A log Ti < 0.02 and 
A log 7b 0.04, and so are on the order of the 
systematics of the analysis (e.g., choice of plasma 
emission model). 

2.2.2. f keV 

A large-scale deabsorption technique was used 
in Snowden et al. (1997) to demonstrate the exis- 
tence of the X-ray emitting bulge with a temper- 
ature of logT = 6.6 which dominates the galactic 
center and the Loop I region, but does not affect 
most of the high galactic latitude sky. 

Pietz et al. (1998) took an interesting, though 
Procrustean, approach to determining the TAE in 
the | keV band. They applied simultaneous spec- 
tral fits to three deep PSPC images from the same 
region of the sky, each with a different absorbing 
column, to determine the temperature of the lo- 
cal and distant components. This method min- 
imizes the problems with inadequate non-cosmic 
background subtraction and avoids the degener- 
acy described in §2.1. 

They then used the fitted temperature for the 
TAE to fit a single-temperature hydrostatic equi- 
librium model of the galactic halo to the observed 
| keV distribution. Given 1) that the high galac- 
tic latitude | keV emission is smooth, 2) that one 
does not expect the canonical model to work well 
near the plane, where the classical “infill” prob- 
lem 5 remains unsolved, 3) that the flattened halo 
model produces very little variation at high galac- 
tic latitude, and 4) that the flattened halo model 
has its greatest variation in the direction of galac- 
tic center where the multiple emission components 
are very confused, it is not surprising that the 

5 The “infill” problem is as follows: at | keV and 1.5 keV 
the average emission level in the galactic plane is roughly 
the same as the average emission at high galactic latitude, 
despite the fact that one would expect galactic disk to 
“shadow” the extragalactic flux. See, for example, Sanders 
et al. (1983). Thus, there must be an additional galactic 
source of diffuse § keV and 1.5 keV emission that almost 
exactly compensates for the emission that is absorbed by 
the galactic disk. However, that source of emission remains 
unkown. 
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Fig. 5. — A demonstration of the bias due the 
method of Snowden et al. (1998). Input tem- 
peratures for both local and distant components 
were used to calculate a model of the observed flux 
as a function of absorbing column density. Those 
model data were then analyzed using the method 
of Snowden et al. (1998) to find “derived” temper- 
atures for the local and distant components. The 
difference between the input and “derived” tem- 
peratures measures the bias in the Snowden et al. 
(1998) analysis. The axes are the input tempera- 
tures, the contours are the “derived” temperature. 
Top: The “derived” temperature of the distant 
component as a function of the true local and dis- 
tant temperatures. The shaded region produces 
measured R2/R1 ratios that do not correspond to 
any temperature. Bottom : The “derived” tem- 
perature of the local component as a function of 
the true local and distant temperatures. Making 
a bad initial assumption will produce an incorrect 
(biased) measure of the distant temperature, but 
iteration will produce the correct distant temeper- 
ature. 


8 



Pietz et al. (1998) model is a “good fit” at § keV. 
Their fitting has found the high galactic latitude 
normalization, but is not good at constraining the 
variation with latitude, since the region with the 
greatest lever arm is also the region that must be 
excluded due to galactic absorption and confusion 
with Loop I. 

Pietz et al. (1998) used a very limited number 
of PSPC pointings towards a fairly atypical re- 
gion (Draco) to make a local/distant separation 
at \ keV. They then claim that the distant com- 
ponent of the \ keV flux (excluding regions with 
High Velocity Clouds or known structures such as 
Loop I) is produced solely by the extrapolation to 
£ keV of the | keV emission produced by their 
galactic halo model. Their evaluation of the de- 
gree to which this extrapolation accounts for the 
angular structure of the - keV emission is overly 
optimistic; as they have no LHB model, they can 
always attribute poorly fitting regions to the ef- 
fects of the LHB. 

2.2.3. Synthesis 

In Snowden et al. (2000), a large (but by no 
means comprehensive) catalogue of 378 ROSAT 
t keV band shadows 6 was used to determine the 
values of the foreground and background emission 
in the R1 and R2 bands for many very localized 
regions at high galactic latitude for both the north- 
ern and southern galactic hemispheres. The | keV 
band emission was determined for the same re- 
gions by subtracting the absorbed EPL from the 
observed emission and deabsorbing the remainder. 
It was found that although both the distant £ keV 
and | keV band emission varied by a factor of five, 
there was no clear correlation between the distant 
j keV emission and the distant | keV emission 
(see the Figure 15 from Snowden et al. 2000). The 
lack of correlation was used to support the argu- 
ment made in Snowden et al. (1998) that the dis- 
tant component that is responsible for the bulk of 
the t keV emission is not the same thermal com- 
ponent that is responsible for the ~ keV emission. 

6 An X-ray shadow is traditionally defined as a region that 
has a deficit of X-rays and the region of the deficit is spa- 
tially matched to an H I 21 cm or 100 ^m emission feature. 
Snowden et al. (2000) expand that definition to regions 

with an excess of X-rays, where the excess is matched to 
the decrease in some measure of the absorbing column. 


2.3. Summary 

Spectral fitting methods are incapable of deter- 
mining the parameters of the thermal components 
due to absorption uncertainty, poor non-cosmic 
background subtraction, and inadequate spectral 
resolution. The two large-scale deabsorption stud- 
ies, executed in different bands, are each blind to 
the analysis in the other. Snowden et al. (1998) 
posited two components to the TAE, but for his- 
torical reasons deferred the necessary analysis in 
the | keV band to expand the proof of their case. 
Pietz et al. (1998) attempted to extrapolate their 
| keV results into the \ keV band, but without 
taking into account the structure of the LHB, or 
the possibility of additional halo components. 

The Snowden et al. (1998) self-consistent so- 
lution for the temperature of the distant compo- 
nent derives a low temperature distant compo- 
nent (logX 1 ~ 6.07) 7 that produces little | keV 
emission. Assuming that both the | keV and the 
| keV distant emission are produced by the same 
thermal component requires a higher temperature 
(log T > 6.14) but, as shown in Figure 4 and § 2.1, 
the higher temperature is inconsistent with the ob- 
served angular variation of the emission with ab- 
sorbing column. Further, the extensive shadowing 
study by Snowden et al. (2000) shows that the 

distant f keV emission is not correlated with the 

f 

distant ^ keV emission. 

Therefore, the distant emission, the TAE, must 
have more than a single thermal component and, 
in fact, may have multiple thermal components. 
The temperature of those components may be spa- 
tially variable. The intensity is undoubtedly vari- 
able. Due to the limitations of the data, no more 
than an incremental improvement can be made by 
extending the Snowden et al. (1998) analysis to 
the R4-R7 bands. We will assume the TAE to be 
composed of two thermal components, a soft ther- 
mal component with temperature Ts dominating 
the ~ keV band, and a hard thermal component 
with temperature T// dominating the | keV band. 
We will derive model temperatures for those com- 
ponents, and will evaluate the extent to which fur- 

7 The temperatures quoted here from Snowden et al. (1998) 
and Snowden et al. (2000) are not quite the temperatures 
found in those papers. The values here have been recalcu- 
lated from the R2/R1 band ratios using an updated detec- 
tor response matrix, while the Snowden et al. papers used 
earlier versions to remain consistent with previous papers. 


9 



ther thermal components may exist. 

3. Method 

Ideally, to determine the thermal and spatial 
structure of the TAE, one would use the angular 
structure of the absorbing column to separate the 
emission in each band into local and distant com- 
ponents, and then one would determine what com- 
bination of thermal components would produce 
the spectrum (or SED) of the distant emission. 
Unfortunately, the separation of the emission of 
any band into local and distant components re- 
quires knowledge of the spectrum of the distant 
component in order to calculate the appropriate 
effective cross-sections. Thus, any determination 
of the nature of the distant emission requires iter- 
ative solution of both the separation into local and 
distant components and the modeling of the ther- 
mal component structure. The extent to which 
such a process can be successful is strongly lim- 
ited by the interplay of the response of the instru- 
ment, the spectrum of the emission, the energy 
dependence of the absorption cross-sections, and 
the available range of absorption column density. 

In the following sections, we will briefly describe 
the available X-ray data and some considerations 
concerning the measurement of the absorbing col- 
umn. Section 3.2 describes the limitations to the 
separation of the X-ray data into local and dis- 
tant components that were alluded to above. Sec- 
tion 3.3 describes the external constraints that can 
be used to set the non-TAE parameters in Equa- 
tion 2. Finally, we detail the iterative method used 
to separate the local and distant components and 
to determine the thermal structure of the distant 
emission. 

3.1. Data 

3. LI. The RASS 

The X-ray data are the single-band RASS maps 
of Snowden et al. (1997), which have 12' x 12' pix- 
els. Of these band maps, the combined R1 and 
R2 bands (the R12 band) map forms a keV” 
map, the combined R4 and R5 bands (the R45 
band) map forms a keV” map, and the com- 
bined R6 and R7 bands (the R67 band) map forms 
a “1.5 keV” map. Histograms of the statistical 
significance of the pixels in each band are given in 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 

tf*.u/Rate 

Fig. 6. — Histograms of the statistical signifi- 
cance of individual 12' x 12' pixels in individual 
and combined bands of the RASS. For bands R4- 
R7, only pixels in the north polar cap (b > 55°, 
45° < £ < 270°) were used. In each panel, the 
solid line is the combined band, the dotted line is 
the lower energy band, and the dashed line is the 
higher energy band. 
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Figure 6. Given the long temporal baseline pos- 
sible with the RASS , the LTE and SSX removal 
were far better than is possible with any single 
pointed PSPC observation 8 . Some systematics re- 
main. Even in seemingly uniform regions, a high 
contrast stretch of the raw RASS images will re- 
veal the scan pattern, more strongly at ~ keV than 
at j keV . 

The R1 and R2 band maps have had point 
sources removed to a level of (R1 + R2) =0.025 
counts s _1 while the R4-R7 band maps have had 
point sources removed to a level of (R4+ R5)=(R6-hR7) 
counts s _1 . Correction of the RASS bands for the 
band- variable point source removal is discussed in 
§3.3.1. 

For reasons that will become apparent shortly, 
the analysis was restricted to 6 > 55°, 45° < I < 
270°; the selection in galactic longitude was set 
to exclude the galactic center/Loop I emission. 
Within this region we also removed pixels con- 
taminated by extended sources, Coma, Virgo, and 
A1367. 

3.1.2. The Absorbing Column 

We have taken as a tracer of the absorbing col- 
umn the IRAS 100 pm maps of Schlegel et al. 
(1997), as no other tracer of the absorbing col- 
umn has either all-sky coverage, or angular reso- 
lution comparable to that of the RASS. As has 
been discussed in Arendt et al. (1998) 100 pm 
emission traces, to varying degrees, the molecular, 
neutral atomic, and ionized atomic components of 
the ISM. (For a review see Chapters 3 k 5 of Kuntz 
2000 .) 

Low velocity disk H I appears to be well traced 
by the 100 f/m emission, but the 100 pm emissivity 
of the H I does appear to vary (Reach et al. 1998). 

The high and intermediate velocity H I appears to 
have very low 100 pm emissivity, and is essentially 
invisible at 100 pm, meaning that the IRAS 100 
pm maps in regions with intermediate (—100 < 
v < —75 km s _1 ) or high (u < —100 km s _1 ) 
velocity clouds will underestimate the absorbing 
column (Wakker &: Boulanger 1986). 

s Removal of residual background contamination for individ- 
ual PSPC pointings after reduction with the ESAS software 
(Snowden & Kuntz 1998) can only be done by reference to 
the RASS . The method is described by Kuntz Snowden 
(1998). 


Only upper limits to the 100 pm emissivity of 
the H II have been determined, so the 100 pm 
emission should be considered a poor tracer of the 
ionized gas (Arendt et al. 1998). 

The emissivity of the molecular gas is a very 
strong function of the cloud temperature. Schlegel 
et al. (1997) attempted to correct for the cloud 
temperature in order to produce maps of the true 
dust column density, but the accuracy of this cor- 
rection, denoted by A', is not well understood. 
Therefore the absorbing column of molecular gas 
0.(}8ven if one assumes a constant dust-to-hydrogen- 
nucleon ratio) may be poorly determined by the 
IRAS 100 pm emission. However, molecular gas 
covers only a small portion of the high galactic 
latitude sky, appearing only where the H I column 
is greater than 3 x 10 20 cm -2 (Reach et al. 1998), 
and the H I 21 cm /IRAS 100 pm ratio is quite 
sensitive to the presence of molecular material. 

The dominant tracer of absorption is H I. (At 
~ keV He provides roughly half of the physical 
absorption, while at energies > 0.5 keV metals 
provide nearly all of the absorption.) The amount 
of absorption due to H II is a matter of debate 
(Arabadjis & Bregman 1999; Laor et al. 1997); 
the cumulative absorbing column due to H II ap- 
pears to be less than 0.2 that of the absorbing col- 
umn due to H I at j keV and may be significantly 
smaller. 

Since the 100 pm emission is very good at trac- 
ing the H I, and the H I 21 cm/100 pm ratio is 
very sensitive to the presence of Ho, the 100 pm 
emission is the single best tracer of the total ab- 
sorbing column. That the 100 pm emission is not 
sensitive to H II is offset by the X-ray absorption 
being less sensitive to H II. To convert from the 
native units of dust emission (MJy sr'" 1 ) to ab- 
sorbing column, we use the ratio of A r (H I)/XVioo 
determined from the northern galactic hemisphere 
(6 > 45°) where ^(H I)< 3 x 10 20 cm -2 , a value 
with varies by ~ 5%, depending upon the latitude 
and column density limits used. 

We have removed the regions where Nh with 
v < — 75 km s -1 has column densities above 1 x 
10 19 cm -2 . This selection removes the strongest 
intermediate and high velocity clouds, which are 
known to be deficient in IRAS 100 pm emission. 
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Fig. 7. — The column densities at various energies necessary to produce a given optical depth. The column 
density marked on each bar is in units of 10 20 cm” 2 . We have assumed a log T — 6.05 thermal background 
in band R12, and a logT = 6.35 thermal background for both R45 and R67. 


3.2. Limits 

Optical Depths: In order to separate the distant 
emission from the local emission, one must have 
absorption over a wide range of optical depths 
within a limited angular region. (The region must 
be small in order for the assumption of a spatially 
invariant local and distant X-ray emission made by 
Equation 1 to be valid.) In practice, the smallest r 
must be < 1, in order to produce a good measure 
of the distant emission, and the largest r must 
be > 2 in order to produce a good measure of //,, 
the local emission. Either of these limits can be 
relaxed if the other end of the optical depth range 
extends well beyond its limit. If the first condition 
is not met, but the latter is, Id may be estimated 
from Id = {Io — Il)c t , where Io is the observed 
flux. 

Since the absoption cross section, a oc E~ 8 ^ 3 , 
regions fulfilling the r requirements at ^ keV gen- 
erally do not fulfill the r requirements at | keV 
or 1.5 keV (see Figure 7). Separation at ~ keV is 
usually possible if Nh < 4 x IO 20 cm -2 , whereas 
separation at | keV is generally only possible if 
Nh » -1 x IO 20 cm” 2 . Thus, while | keV sepa- 


ration into Id and //, is usually possible at high 
Galactic latitude, | keV separation is generally 
only possible at low latitude. 

Count Rates: The difficulty in separating Id 
from I 1 at | keV is exacerbated by the | keV 
count rates, which are typically one-tenth of the 
count rates at | keV (note the lower statistical sig- 
nificance of higher energy band pixels in Figure 6). 
Thus, separation of the RASS is generally possi- 
ble for the R1 and R2 bands, individually, but at 
higher energy, it is only possible in the combined 
R45 and R67 bands. 

Regions with sufficiently high column densi- 
ties to show good | keV shadows are typically 
near the galactic plane, where there are likely 
to be many emission components (and absorbing 
screens) along a given line of sight. Thus, although 
there are regions for which the local/distant sep- 
aration is mathematically possible, those regions 
are unlikely to produce a meaningful measure of 
either the local or the distant emission. 
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3.3. External Constraints 

3.3 A. Properties of the EPL 

The EPL is the aggregate of unresolved sources. 
It must be the result of a complex convolution 
of source populations with different spectra, the 
distribution of redshifts, and the distribution of 
intrinsic absorption. None of the existing popu- 
lations of point sources thought to compose the 
EPL has spectra as hard as that of the EPL. 
Given the cosmological significance of the EPL, a 
number of groups have attempted to measure the 
photon index (T) and normalization (At) of the 
EPL (i.e., the spectrum of the aggregate of unre- 
solved sources). Indices measured from ROSAT 
data alone have been higher than those measured 
from ASCA , which provides a significantly longer, 
absorption free, lever arm (e.g. Chen et al. 1997). 
It has been suggested by a number of authors that 
this discrepancy is due to calibration uncertain- 
ties in the ROSAT PSPC. Fits using both ROSAT 
and ASCA data tend to indices more similar to 
the ASCA values. We have chosen the values of 
Chen et al. (1997), T = 1.46, N T = 10.5 keV 
cm" 2 s -1 sr -1 keV -1 , as being representative of a 
number of recent measures of the EPL from joint 
ROSAT/ ASCA fits. 

The Chen et al. (1997) point source removal 
limit for the region from which the T = 1.46, 
Nr = 10.5 EPL was derived was ~ 0.00524 counts 
s -1 for the R47 band. To correct to the point 
source detection limit of the RASS , nearly an or- 
der of magnitude larger, we used the log N — log 5 
relation of Hasinger et al. (1998) and have assumed 
a photon index for the point sources between the 
point source removal limits of Chen et al. (1997) 
and the RASS . We have assumed that those point 
sources can be characterized by the point sources 
detected by Chen et al. (1997) at somewhat lower 
fluxes (< 0.00524 counts s -1 but > 0.0003 counts 
s -1 ). A T = 2.54 index was derived by Chen et al. 
(1997) for a sample of confirmed QSO’s detected 
in this flux interval. However, the confirmed QSO 
sample is softer than the complete sample of ob- 
jects in this flux interval, so a harder index may be 
more apropriate. We have performed the analysis 
using both T = 2.54 and V = 1.46 for the point 
source correction 9 . 


9 The point source correction for T = 2.54 requires a nor- 


Caveat: Although much effort has been ex- 
pended to determine the spectra of individual 
sources, the spectrum of the unresolved extra- 
galactic background, the EPL, is unknown in the 
1 keV band. We have assumed that the power 
law that describes the extragalactic background 
for E > 1.0 keV extends into the \ keV band. 

It has been suggested that the spectrum of the 
extragalactic background turns up below 1 keV. 
However, that turn-up can not be too great; Snow- 
den et al. (2000) show that the isotropic extra- 
galactic background at \ keV cannot be more than 
~ 325 x 10” 6 counts s -1 arcmin -2 . If the spec- 
trum breaks at 1 keV, then T < 1.64; if it breaks 
at 0.75 keV, then T < 1.67. Thus, even if there is 
a turn-up in the extragalactic spectrum, the effect 
is small. 

3.3.2. | keV Emission from the LHB 

Given all of the caveats placed on | keV lo- 
cal/distant separation in §3.2, any external con- 
straint on the magnitude of the local | keV emis- 
sion is useful. Lines of sight useful for measur- 
ing the local component of the | keV emission 
must be chosen carefully. The highest signal-to- 
noise measure of the local | keV emission remains 
that of Snowden et al. (1993) towards MBM 12, 
(a molecular cloud situated just inside the LHB) 
which yields 2 a upper limits of 23 x 10 -6 counts 
s -1 arcmin -2 (| keV) and 39 x 10 -6 counts s -1 
arcmin"* 2 (1.5 keV). In this direction, the - keV 
flux is 385 x 10~ 6 counts s -1 arcmin"" 2 . 

If the temperature of the LHB is ~ log7T = 
6.10, then the R45/R12 ratio is 0.017, and the 
| keV count rate due to the LHB in the direction 
of MBM 12 is 7.1 x 10 -6 counts s -1 arcmin -2 , well 
below the upper limit and similar to the measured 
value of ~ 5 x 10 -6 counts s -1 arcmin -2 . At high 
galactic latitudes, where the LHB flux is on the 
order of 400-600 x 10 -6 counts s -1 arcmin -2 , the 
| keV counts due to the LHB will be 8.75 — 13.25 x 
10 -6 counts s -1 arcmin -2 , which is a substantial 

malization of 8.032 x 10“ 13 erg cm -2 s -1 deg -1 in band 
R47 and produces the corrections Rl = — 45.9, R2=— 38.7, 
R4=4.3, R5=5.1, R6=6.6, R7=3.2, where the values are 
in units of 10 -6 counts s _1 arcmin -2 to be added to the 
EPL. The point source correction for F = 1.46 requires 
a normalization of 9.712 x 10~ 13 erg cm -2 s -1 deg -1 in 
band R47 and produces Rl = — 14.3, R2= — 15.3, R4=3.2, 
R5=4.9, R6=8.4, R7=5.9. 
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fraction of the high galactic latitude | keV flux, 
typically ~ 35 x 10~ 6 counts s~* arcmin -2 . Thus, 
ignoring the LHB contribution at - keV in the 
manner of Pietz et al. (1998) is unwise. 

3.4. The Iterative Algorithm 

We use a bicyclic iterative scheme to solve for 
the distant component in each band and, simul- 
taneously, fit two thermal components to the re- 
sulting SED. The process is described graphically 
in Figure 8. Given the difficulty of performing 
a local/distant separation in the R4-R7 bands, 
we have treated those bands somewhat differently 
than the R1 and R2 bands; the Rl and R2 bands 
are separated into local and distant components 
with a method similar to that used by Snowden 
et al. (1998), while the R4-R7 bands are corrected 
only for galactic absorption after the contribution 
of the LHB has been removed. The lack of cor- 
relation between the | keV and | keV bands, as 
well as the low model temperature found for the 
^ keV emission by Snowden et al. (1998) suggest 
that there is little contribution from the \ keV 
emission to the | keV emission (when compared 
to the much greater flux at ~ keV ) and vice versa. 
The fitting algorithm begins with this assumption, 
but the assumption is not reified by the fitting al- 
gorithm. 

The first cycle (the left side of Figure 8) iter- 
atively determines the temperature of the softer 
component from the bands Rl and R2, after sep- 
arating the local from the distant emission. The 
second cycle (the right side of Figure 8) iteratively 
determines the temperature of the harder compo- 
nent from the bands R4 through R7. Each cycle 
passes to the other the extrapolation into the com- 
plementary band range of the component it has fit. 
The first cycle finds the single best thermal com- 
ponent fit to the Rl and R2 bands, and extrapo- 
lates that component into the R4-R7 bands. That 
extrapolation is then subtracted from the R4-R7 
data, and a single thermal component is fit to the 
remainder. That thermal component is then ex- 
trapolated to the R1-R2 band range, where it is 
subtracted from the R1-R2 data. The process is 
then repeated. 


3.4.1. Details of Local/Distant Separation 

The process by which the distant emission is 
separated from the local emission in bands Rl and 
R2 is similar to that used by Snowden et al. (1998). 
Both the local and distant emissions are assumed 
to be spatially invariant over 4? 8 x 4? 8 regions. 
The bulk of the distant emission in the R12 band 
is assumed to come from a single thermal compo- 
nent. From each /MSS single band image is sub- 
tracted the contribution of the EPL, after absorp- 
tion by the galactic ISM, and the contribution to 
that band from the hard thermal component, also 
after absorption by the galactic ISM. The X-ray 
data from each 4? 8 x 4? 8 region are then binned by 
the absorbing column density into 1 x 10 19 cm -2 
wide bins. The binned data are then fit by Equa- 
tion 2, assuming some logT for the distant emis- 
sion. The Lampton et al. (1976) criteria are used 
to determine the uncertainties in the fit parame- 
ters. The fitting is done for each 4? 8 x 4? 8 region 
for which the Nh < 4 x 10 20 cm" 2 . The 4?8 x 4?8 
regions overlap by 1?6, so the solutions are applied 
only to the center ninth of the region (1?6 2 ). 

In their work, Snowden et al. (1998) found that 
the uncertainty in the distant component could be 
substantially reduced by taking the quantity 

J^ = (/o-7i> r (4) 

rather than the value produced by the fit. We have 
followed this procedure, although, owing to slight 
differences in the fitting procedure, we find that 
the improvement is less significant. 

Once the separation is accomplished for both 
the Rl and the R2 band images, the hardness 
ratio is determined for the distant component. 
This leads to a new estimation of the logT? for 
the distant emission, and the process is repeated 
until the assumed temperature and the temper- 
ature recovered from the hardness ratio agree to 
A log Ts < 0.005. Convergence typically occurs in 
three iterations. 

3. 4 .2. Details of the R4-R7 SED Fitting 

From each RASS single band image is sub- 
tracted the contribution from the EPL, absorbed 
by the galactic ISM, the contribution from the soft 
thermal component that dominates the Rl and R2 
bands, absorbed by the galactic ISM, and the con- 
tribution by the LHB, without absorption. The re- 
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Pig g — Flowchart for the bicyclic iterative fitting of two thermal components to the SED of the TAE. 
T s is the model temperature of the softer component (which is fit in the left iterative cycle) and Th is 
the model temperature of the harder component (which is fit in the right iterative cycle). R n is an image 
of the observed emission in band n, E n is an image of the emission extrapolated from other energies into 
band n, I n is an image of the TAE emission in band n, corrected for absorption, P n is the value of the 
extragalactic power law in band n, and C is the absorbing column. The subscripts L and D denote the local 
and distant components. B n is the band value derived from a Raymond & Smith spectrum. The U T - W 
is a consistency check for the temperature of the component. 


mainder is smoothed by a three pixel boxcar and 
is then deabsorbed (i.e., multiplied by e T where r 
is the effective cross-section multiplied by the col- 
umn density) assuming some log Th for the distant 
emission. 

After all four bands have been deabsorbed in 
this manner, the mean deabsorbed SED for a re- 
gion with b > 55° and 45° < l < 270° is calcu- 
lated. A thermal model is fitted to this SED to 
determine the new log T//, and the process is re- 
peated until A log T h < 0.005. Convergence typ- 
ically occurs in three iterations. Convergence be- 
tween the two halves of the iteration usually occurs 
in three or four iterations. 


4. Analysis of Current Results 

The results of the local/distant separation and 
the SED fitting are shown in Figure 9. The top 
images show the soft and hard components in the 
bands in which they dominate, and below that 
are their uncertainties (3 a for the soft, la for the 
hard). The third row contains the spatial dis- 
tribution of the temperature of the two compo- 
nents. For the soft band the temperature is in 
the form of the R2/R1 band ratio, while for the 
hard component the temperature is the fitted Ray- 
mond & Smith (1977) equilibrium temperature. 
The fourth row contains the absorbing column, 
and the raw R47 band image. 
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Fig. 9. — Previous page: The spatial distribution of measured quantities. All images are of the northern 
galactic hemisphere. Galactic longitude 0° is down, 90° is on the left. The longitude spacing is 30°, the 
latitude spacing is 15°. Note that the pixelization is 1?6. a: The soft TAE component as seen in the R12 
band. Greyscale: minimum 0, maximum 5 x 10“ 3 counts s" 1 arcmin -2 . b: The hard TAE component as seen 
in the R47 band. Greyscale: minimum 0, maximum 5 x 10 -4 counts s"" 1 arcmin~ 2 . c: The 3a uncertainty of 
the soft TAE component. Greyscale as in a. d: The lcr uncertainty of the hard TAE component. Greyscale 
as in b. e: The R2/R1 band ratio for the soft component. Greyscale: minimum 0, maximum 3. f: The 
log of the fitted temperature of the hard TAE component. Greyscale: minimum 6.0, maximum 7.0. g: The 
absorbing column density. Greyscale: minimum 0. maximum 4 x 10 20 cm -2 , h: The observed R47 band 
flux. Greyscale as in b. Note that the bulk of the observed structure is due to the scan paths. 


4.1. The TAE 

4-1- 1- The Soft Component 

The temperature of the soft component was de- 
rived from the R2/R1 ratio over the north polar 
cap region ( b > 45° , 45° < t < 270°, exclud- 
ing regions with strong X-ray clusters of galaxies 
or strong (TV// > 10 19 cm -2 ) high and intermedi- 
ate velocity clouds). The histogram of the R2/R1 
band ratio over the north polar cap region is shown 
in Figure 10, along with the conversion to Ray- 
mond & Smith model temperatures. The FWHM 
of the R2/R1 ratio is 0.670-1.190, and the distribu- 
tion is asymmetric, with a strong high R2/R1 tail. 
The values in the high tail do not correspond to a 
Raymond & Smith model temperature, and thus 
indicate that the bulk of the dispersion is due not 
to an intrinsic spread of temperatures, but to the 
uncertainty of the local/distant split. The mean 
uncertainty in R2/R1, obtained by propagating 
the uncertainties in the distant components of R2 
and R1 is sufficient to produce the entire observed 
FWHM. Unlike Snowden et al. (1998), we have 
not smoothed the distant emission images before 
calculating the R2/R1 ratios, so our uncertainties 
are somewhat higher. 

In order to minimize the effects of the tail we 
have chosen to characterize the R2/R1 ratio by 
the mode. The mode is 0.950; the correspond- 
ing model temperature is logXs = 6.056. Snow- 
den et al. (1998) found logT = 6.09 (FWHM = 
6.01 - 6.17) for the soft distant component, with- 
out an extrapolation of a hard thermal compo- 
nent. If we follow Snowden et al. (1998), and ig- 
nore the extrapolation of the hard component, we 
find an R2/R1 ratio of 0.92, with a correspond- 
ing model temperature of log T = 6.040. Since we 
used a much more restricted portion of the sky, 


and a different fitting algorithm, the difference be- 
tween the value from Snowden et al. (1998) and 
our value (without the extrapolation of the hard 
component) is a good measure of the systematic 
uncertainties. 

Systematics: Most of the systematics in the 

local/distant separation have been discussed in 
Snowden et al. (1998). The magnitude of the dis- 
tant component is not correlated with the absorb- 
ing column, though the uncertainty will be. The 
value of the R2/R1 ratio for the distant compo- 
nent is not correlated with the value of R2, Rl, or 
the magnitude of the absorbing column. Snowden 
et al. (2000) discuss the variation of the R2/R1 ra- 
tio of the distant component with galactic latitude 
and longitude (they find none), and although we 
have a greater number of pixels, Snowden et al. 
(2000) have the advantage because their regions 
were larger, and were carefully chosen to have a 
large dynamic range in absorption 

Although the R2/R1 ratio of the distant compo- 
nent is not correlated with A r //, the derived tem- 
perature of the distant component may still be 
systematically affected by the absorbing column. 
If the absorbing column is systematically underes- 
timated, one might overestimate the model tem- 
perature, and that overestimate may depend upon 
the absorbing column. As can be seen in Fig- 
ure 11, for the approximate background temper- 
ature and range of columns used ( Nh < 4.0 x 
10 20 cm" 2 ), the overestimate has little depen- 
dence on absorbing column. If Nhii/Nhi ~ 0.33 
(Reynolds 1991) and (jhii/ohi ~ 0.62 (assum- 
ing that N Hell /A//e ~ 0.5 as per Domgorgen 
& Mathis 1994) then the absorbing column due 
to H II will be ~ 0.2 times that due to H I. 
Assuming a distant component temperature of 
logT — 6.15, if one does not account for the 20% 
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extra absorption due to the H II, one would over- 
estimate the distant component temperature by 
A log 0.05. 

One might also expect to find that the lack 
of a sufficient dynamic range in the absorption 



R2/R1 


Fig. 10. — The distribution of the R2/R1 band 
ratios for the LHB and TAE. The squiggle above 
the histogram shows the conversion to Raymond &; 
Smith model temperature. For the TAE, both the 
R2/R1 ratio and Raymond & Smith conversion are 
after deabsorption. For the LHB the histogram is 
of the R2/R1 ratio derived from I ( L values; the 
temperature conversion assumes a 4 x 10 18 cm" 2 
absorption. Note that only the lower branch of the 
temperature conversion is of interest as the upper 
branch would over-predict the R47 flux. 


might produce a systematically high (or low) esti- 
mate of the distant component. We recalculated 
the R2/R1 ratio using only those 1?6 2 pixels with 
“good” fits, defining “good” fits to be those for 
which the ler lower limit was greater than zero. 
The difference between this value of R2/R1 and 
that calculated from all of the pixels was insignif- 
icant. 

4-1.2. The Hard Component 

The temperature of the hard component was 
derived from fitting the R4-R7 SED. The uncer- 
tainty of the value was determined from the dis- 
persion of the temperatures found from fitting the 
SED of each 1?6 2 region. The histogram of the dis- 
tribution of temperatures is shown in Figure 12. 
The mode of the temperature distribution was 
log Tn = 6.42 and the FWHM of the distribution 
was 6.30-6.56. 

To within the uncertainties, the spatial distri- 
bution of the temperature of the hard component 
is uniform. The gradient from bottom to top in 
Figure 9 matches the scan pattern seen in the 
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Fig. 11. — The derived R2/R1 ratio as a function 
of Nh if Nh were underestimated. We assumed 
the distant emission to have logT = 6.15 and mea- 
sured the R2/R1 ratio of the deabsorbed emission 
when the absorbing column had been underesti- 
mated by (from bottom to top) 0% to 30% in steps 
of 5%. 
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RASS images; the pattern has an amplitude of 
A logT — 0.1. 

4-1.3. Comparisons 

Comparison of the components: If the TAE 
were produced by a single component, then one 
would expect the R12 and R45 bands for the TAE 
to be correlated. The lack of correlation is shown 
in Figure 13. 

Comparison of one- and two-component mod- 
els: Given the R1 and R2 measures of the TAE 
(from the local/distant separation), and the R4- 
R7 measures of the TAE (from direct RASS ob- 
servation), we can construct the SED for each 1?6 2 
region on the sky. We can then measure the rel- 
ative goodness of fit for both the one- and two- 
component models. 

To minimize the introduction of a bias through 
deabsorption, we created an absorbed SED from 
that part of the sky with 0.95 < A r ///10 2 ° < 1.15 
cm -2 . For all bands we have taken the SED ele- 
ments to be Io — ( II + Ie) where Ie is the contri- 
bution due to the absorbed EPL. We have then fit 
the two component model to each 1?6 2 pixel using 
the previously derived temperatures, and deter- 
mined the mean \ 2 per fit, \l = 3.39 (v = 4). 

The temperature of the single-component model 
was determined iteratively in a manner similar to 
that used for the two-component model. We fit a 
single-component model to the SED used for the 
two component fit. That temperature was then 
used to do the R1 and R2 band local/distant sep- 
aration, and to redetermine the SED. The single- 
component model was refit to the R1-R7 SED, 
and the procedure was repeated until sucessive 
fits produced the same value. We then found 
the temperature for which a single-component 
model had the best mean x 2 per fit, x 2 = 6-62 
( 1 / = 5) at logT = 6.20. If the single-component 
model were correct, one would expect that the 
temperature produced by fitting the SED should 
be consistent with that found from the R2/R1 
ratio. Using logT = 6.20 as the TAE tempera- 
ture in the local/distant separation, we find the 
unabsorbed R2/R1 band ratio to be 0.97, equiv- 
alent to a Raymond & Smith model temperature 
of logT = 6.07. Therefore, although our proce- 
dure attempts to make the temperatures found by 
SED fitting and by local/distant separation con- 



Fig. 12. — The distribution of the temperatures 
produced by fitting Raymond & Smith models to 
the R4-R7 SED after the emission due to the LHB 
and the soft component of the TAE have been re- 
moved. 
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Fig. 13. — Correlation of the | keV (R12) band 
emission of the TAE and the ~ keV (R45) band 
emission of the TAE. The points are shown with- 
out error bars for clarity. 
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Fig. 14. — Fitting one- and two-component TAE 
models to the R1-R7 band SED. Both panels use 
a “mean” SED calculated from all regions where 
0.95 < Nff/ 10 20 < 1.15 and (45° < i < 270°, 5 > 
55°). The uncertainties are representative for in- 
dividual 1?6 2 pixels. Top: The data are the x, 
the solid line is the best fit two-component TAE 
model, the dotted line is the hard component, the 
dashed line is the soft component. Bottom: The 
data are the x, the solid line is the best fit one- 
component TAE model. Note that the SED used 
for the one-component model is slightly different 
from that used for the two-component model. See 
text for explanation. Note that the relative x 2 val- 
ues given in the text were derived from all of the 
pixels within (45° < t < 270 °, b > 55°) while the 
fits shown are to the mean of a subsample of the 
pixels used. 


sistent, no single-component model can be made 
consistent. 

The two-component TAE model is clearly a sig- 
nificantly better fit. Results for a “mean” SED 
derived from a limited range of absorbing column 
density (0.95 < A 7 //IO 20 < 1.15) can be seen for 
both models in Figure 14. The single-component 
temperature derived here is consistent with other 
single-component analyses. The relatively poor fit 
of the single-component model is due to forcing 
the SED fit to be consistent with the local/distant 
separation. 

More Components?: The data are clearly insuf- 
ficient to determine the existence of further com- 
ponents; any combination of thermal components 
that mimic the two derived here will work as well. 
The spatial distribution of the temperature of the 
hard component (Figure 9) suggests that that tem- 
perature may vary with position. Whether there 
is a single component with a spatially variable 
temperature, or several components with different 
temperatures must, at this point, remain a seman- 
tic distinction. 

Conversely, one can ask if the hard thermal 
component is actually due to some component of 
the extragalactic background. If so, then the spec- 
trum of the background must turn up at ~ 1 keV, 
but turn back down before the | keV band; pos- 
sible, but contrived. 

It is likely that some portion of the hard ther- 
mal component is due to unresolved Galactic stars. 
Stellar X-ray emission is usually modeled by ther- 
mal plasmas with logT ~ 6.3 - 6.5. A calculation 
of the X-ray flux due to stars, similar to that made 
by Guillout et al. (1996), shows that no more than 
~ 30% of the hard component flux can be due to 
unresolved stellar point sources. At high galac- 
tic latitudes, in directions of nominal diffuse back- 
ground levels, only ~ 25% of the | keV and 1.5 
keV band flux which is not due to the LHB or the 
EPL can be due to unresolved stars. 

Comparison of point source spectrum inodels: 
As noted in §3.3.1, one must correct the RASS 
for point sources below the RASS removal limit 
but brighter than those discarded by Chen et al. 
(1997) when measuring the spectrum of the EPL. 
The spectrum of those sources is poorly under- 
stood in the £ keV and | keV bands. The re- 
sults presented above use T = 2.54. for the point 
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source correction. If a T = 1.46 spectrum were 
used instead, the changes to the derived temper- 
atures is minimal; log Ti — 6.10, log Ts = 6.08, 
and log T h = 6.39 for T = 1.46 vs. logT/, = 6.10, 
log Ts = 6.06, and log I* = 6.42 for T = 2.54. 

Comparison of plasma models: From Table 1 
it will be remembered that the model tempera- 
ture derived depended, to some degree, upon the 
plasma model used. To some extent, this disagree- 
ment is immaterial; no plasma model has been 
able to reproduce the line spectrum of the diffuse 
X-ray background (Sanders et al. 1998; Vallerga 
& Slavin 1998; Deiker et ah 1997). Therefore, the 
model temperatures derived here may be nothing 
more than convenient markers on an artificial but 
universally accessible scale. However, these mark- 
ers certainly allow comparison of some measure of 
the relative temperatures of the gas. 

Using XSPEC (Arnaud 1996) to fit Mewe- 
Kaastra-Leidahl models to unabsorbed Raymond 
& Smith models convolved with the ROSAT re- 
sponse, we find that 

log Tmkl ~ 0.92 -1- 0.83 log Trs • (5) 

Thus, the equivalent MKL temperatures will be 
systematically lower than the Raymond Sz Smith 
values quoted here. 

4.1.4. The Diffuse Background Budget 

Table 2 presents the budget for the diffuse soft 
X-ray background toward the north galactic pole 
( b > 85°, 45° < £ < 270°, with the region around 
the Coma cluster excluded). The line marked ‘'Re- 
mainder” is the total observed flux from which has 
been subtracted the contributions from the LHB, 
the EPL, and the point source correction necessary 
for point sources not excluded from the RASS , but 
too bright to have been included in the Chen et al. 
(1997) formulation of the EPL. Hard and soft ther- 
mal components with log Th = 6.42 and log Ts = 
6.06 were fit to the remainder to determine their 
normalizations. The final line shows what portion 
of the remainder is due to unresolved stellar X-ray 
sources. It should also be noted from the work of 
Snowden et al. (2000) that of the remainder no 
more than 80xl0” 6 counts s _1 arcmin~ 2 in the 
~ keV band can be truly extragalactic. (That is, 
after the removal of the EPL, the minimum flux 
observed originating from beyond the galactic H I 
layer is 80xl0 -6 counts s“ x arcmin” 2 .) 


4-1.5 . Spectral Fitting Redux 

Using the knowledge gained above, one can re- 
turn to the spectral fitting method. We have 
used three deep pointed PSPC exposures to test 
the consistency of our broad-band RASS results 
with more detailed spectral fitting. Whereas Pietz 
et al. (1998) selected pointings that were close to- 
gether but had different absorbing columns, we 
used maps of the distant £ keV emission to choose 
directions with different TAE intensities. We then 
searched the ROSAT archive at the HEASARC for 
pointed observations longer than 10 ks in those 
directions, and which did not include clusters of 
galaxies or particularly bright point sources. Ta- 
ble 3 details these observations. The three obser- 
vations, one exposure of the Lockman Hole deep 
survey, one of an AGN, and one of a star, were 
then prepared for spectral fitting. We first exam- 
ined the light curves of the observations and elim- 
inated times of bright scattered solar X-ray con- 
tamination, short-term enhancements (most likely 
a background due to auroral X-rays), and long- 
term enhancements (the observations were rela- 
tively unaffected by the last). We examined im- 
ages of the data and created an extraction region 
which included most of the area inside of the ring 
of the window-support structure but excluded the 
targets of the AGN and stellar observations. The 
same extraction region was used on all observation 
data sets. We binned the data into pulse-height 
channels with a minimum of 25 counts and cre- 
ated ancillary region files (ARFs). 

XSPEC V.10 was used to fit the three observa- 
tions simultaneously with several different models. 
We first fit the data with the two-component TAE 
model: an unabsorbed thermal spectrum for the 
LHB, two absorbed thermal spectra for the TAE, 
and the absorbed EPL. For the fit we fixed the 
power-law slope to the Chen et al. (1997) value of 
1.46 but allowed the normalization to vary. (The 
fitted normalization was 12.1 keY cm" 2 s'" 1 sr -1 
keV -1 .) The temperatures of the three thermal 
components were allowed to vary but they were 
constrained to be the same for all three data sets. 
The normalizations of the local component and the 
soft absorbed component were fit independently 
for the three data sets but the normalizations of 
the hard thermal component and power law were 
constrained to be the same. The data were fit in 
the 0.12-1.90 keV band. Table 4 lists the results, 
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Table 2 

Diffuse Background Budget 1 


Component 

Observed Fluxes 

Fluxes Before Absorption 5 


10 6 counts s 1 

arcmin - ^ 

10 6 counts s 1 

arcmin 2 


R12 

R45 

R67 

R12 

R45 

R67 

Total Observed a 

1036.6 

122.3 

122.6 




LHB 

575.2 

9.6 

0.4 




EPL 

94.0 

57.1 

105.4 

245.6 

60.5 

107.4 

PSC C 

-31.5 

8.9 

9.6 

-84.6 

9.4 

9.83 

Remainder 

399.0 

46.7 

7.2 




Best Fit Soft d 

349.8 

8.9 

0.3 

1015.6 

10.5 

0.4 

Best Fit Hard e 

50.7 

37.4 

10.1 

124.9 

40.0 

10.5 

Stars f 

9.0 

9.7 

4.3 





A RASS data taken from b > 85°, 45° < l < 270°. 

b Mean absorbing column was 1.32 x lO 20 cm“ 2 . 

c The point source correction described in the text, §3.3.1 

d The best fit emmission measure was 7.3 x 10" 3 cm" 6 pc. 

e The best fit emmission measure was 2.3 x 10" 3 cm" 6 pc. 

f Unresolved stellar X-ray sources which contribute to the hard thermal 
component. 


Table 3 

Pointed Observations Used 


Observation 

N h 

l 

b 

Exposure(ks) 
Total Useful 

rp900029a04 

0.56 x 10 20 

149.51 

53.16 

46.7 

46.7 

rp700996n00 

2.16 x 10 20 

217.00 

55.45 

19.6 

14.4 

rp200721n00 

1.88 x 10 20 

86.55 

40.91 

46.9 

28.0 
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Fig. 15. — Simultaneous fits of three PSPC spectra from deep pointings. Top: Fit by a model containing the 
EPL, LHB and a two-component TAE. 2nd Panel: The residuals to the two-component fit. 3rd Panel: Fit 
by a model containing the EPL, LHB, and a one-component TAE; all the spectra are fitted with the same 
temperature TAE. J^th Panel: The residuals to the one component fit. 5th Panel: Fit by a model containing 
the EPL, LHB, and a one-component TAE; each spectrum is fitted with its own TAE temperature. 6th 
Panel: The residuals to the one-component fit. Bo^om: The best fit two-component models minus the 
best fit one-component models, where the one-component models were required to have the same TAE 
temperature, i.e. the models in the third panel. 



Table 4 
Spectral Fits 


Model logJ L 

A'z. a 

log Tp\ 

AW 1 

logT D2 A’o-/ xl 

V 

2 Comp 6.10 ±0.02 

1.28 ±0.53 

6.00 ± 0.02 

< 0.77 

6.32 ±0.02 0.64 ±0.24 1.22 

450 


0.73 ±0.14 


1.72 ±0.60 




2.61 ±0.04 


4.14 ±0.89 



1 Comp b 6.17 ± 0.03 

1.76 ±0.04 

6.56 ±0.07 

< 0.05 

1.39 

452 


1.28 ±0.04 


0.20 ± 0.09 




2.61 ±0.04 


< 0.05 



1 Comp c 6.10 ±0.02 

1.28 ± 0.24 

6.31 ±0.23 

0.64 ±0.35 

1.26 

450 


0.09 ± 0.14 

6.24 ± 0.02 

1.25 ±0.42 




1.58 ±0.26 

6.13 ±0.01 

3.01 ±0.71 




a Units are 10 3 photons cm 2 s 1 keV 1 in a solid angle of 8.54 x 10 5 sr. 

b The temperature of the absorbed thermal component was constrained to be the same for all data sets. 
The normalization of the extragalactic power law was fixed at the fitted value from the two-component 
model. 

c The temperature of the absorbed thermal component was fitted independently for all the data sets. 


which show a reasonable goodness of fit (\l = 1.22 
with v — 450). The fitted values of the tempera- 
tures are consistent with our broad-band analysis. 

We also fit the data with the one-component 
TAE model: an unabsorbed thermal spectrum for 
the LHB, one absorbed thermal component for the 
TAE, and the absorbed EPL. The temperatures of 
the thermal components and the normalization of 
the power law were constrained to be the same for 
all data sets but the normalizations for the ther- 
mal components were fit individually. With two 
fewer fitted parameters the reduced \ 2 value in- 
creased from 1.22 to 1.38. This rather significant 
reduction in the quality of the fit relative to the 
two-component model supports our claim that the 
one-component model cannot, in general, be fit 
simultaneously to multiple directions on the sky. 
Examination of the fitted parameters (Table 4) re- 
veals a much higher temperature for the absorbed 
thermal component (logT — 6.56), and lower val- 
ues for the normalizations. The fitting algorithm 
forced the local emission to accept a portion of the 
distant emission. 

To examine the situation further, we allowed 
the temperature of the absorbed thermal compo- 


nent in the one-component model to be fit indi- 
vidually to each data set. This produced a consid- 
erable improvement in the \ 2 value, to 1.26, now 
with the same number of fitted parameters as the 
two-component model. However, the fitted values 
for the temperatures (Table 4) were completely in- 
consistent with a constant value (as expected with 
the large change in the quality of the fit). The 
brighter the \ keV TAE the cooler the fitted tem- 
perature. This is exactly the expected result if 
the ~ keV intensity arises from a different emis- 
sion component than the | keV intensity. Thus 
the two-component model is completely consistent 
with the detailed spectral fits of the three obser- 
vations while the one-component model is not. 

The results listed in Table 4 also demonstrate 
the ease with which the fitting algorithm can shift 
intensities and spectral parameters between one 
component and another. The choice of details 
used for the spectral model strongly affects the re- 
sults which are obtained. The two separate fits to 
the one-component model provide an excellent ex- 
ample. In the first case, when the temperature of 
the TAE is constrained to be the same for all three 
data sets, the fitted temperature is very high and 
the fitted intensities are very low; such low intensi- 
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ties are inconsistent with the variation of the halo 
emission which has been determined by shadowing 
experiments. This is clearly the case with the ob- 
servation in the high-intensitv direction (towards 
the Draco enhancement) where the fit provides a 
low value for the upper limit; flux from the TAE 
has been shunted to the LHB. 

The difference between the temperature derived 
for the hot component by spectral fitting and that 
derived by SED fitting should not be considered 
too seriously; the two methods have very different 
systematics. Whereas, in the bicyclic iteration, 
the normalization and temperature of the LHB 
was very strongly constrained by the anticorrela- 
tion with the absorbing column, in the spectral 
fitting, the properties of the LHB are less strongly 
constrained, and so the fitting routine can dis- 
tribute flux among the various components more 
freely, producing the best fit for those three par- 
ticular spectra. 

4.2. The LHB 

Overview: The local/distant split that produces 
the R1 and R2 measures of the TAE also pro- 
vides a measure of the emission due to the LHB. 
We find an LHB temperature of log!/, = 6.10 
(FWHM = 6.02-6.26) from the R2/R1 band ratio 
1.04 (FWHM = 0.890-1.220). Although Snowden 
et al. (1998) and Snowden et al. (2000) contain 
detailed discussions of a similar analysis, as ap- 
plied to the LHB, several points are worth noting. 
First, the temperature of the LHB, as derived by 
this analysis or by spectral fitting, depends sensi- 
tively upon the temperature found for the TAE. 
The total R1 and R2 band emission must be con- 
served, and the imposition of a given TAE temper- 
ature will determine the temperature to be found 
for the LHB. However, there are some checks upon 
the temperature of the LHB. Since, in the galac- 
tic plane, the column density in many directions is 
sufficient to block background ~ keV emission, one 
can determine the R2/R1 ratio for the LHB with- 
out TAE contamination. After marshaling some 
essential data about the LHB, we will compare 
our measure of the LHB temperature to those in 
the literature. 

Shape and Size: The shape of the LHB was de- 
termined by Snowden et al. (1998) from the all-sky 
deabsorption, but the size can not be set without 
knowing the emission measure within the LHB. 


The standard scaling of the LHB is given by the 
observation of the molecular cloud MBM 12, which 
produces a barely significant \ keV shadow, im- 
plying that it is just barely within the LHB emis- 
sion. The distance to MBM 12 has recently be 
redetermined to be ~ 90 pc by Sfeir et al. (1999). 
Thus, the high scaling of Snowden et al. (1998), 
0.23 pc (10 -6 counts s -1 arcmin -2 ) -1 , should be 
used. 

Internal Absorption: Sfeir et al. (1999) has 
compiled a homogeneous data set of 456 lines of 
sight with Na I D absorption spectra with which 
they have mapped the Local Cavity. Their com- 
parison of the Snowden et al. (1998) LHB con- 
tours to the Local Cavity shows that the bulk of 
the LHB fits inside the distance at which the cu- 
mulative column density is significantly less than 
2 x 10 19 cm -2 (Jl\ = 10mA); only some of the 
more peripheral regions of the LHB overlap direc- 
tions/distances with such high column densities. 

Absorption of the LHB emission is produced 
by the Local Interstellar Cloud (LIC) and several 
similar clouds within ~ 5 pc. A mean absorb- 
ing column due to these features can be derived 
from Hutchinson et al. (1998) who showed that the 
mean absorbing column within ~ 50 pc was about 
4 x 10 18 cm -2 , about a factor of two greater than 
the highest LIC column densities (Lallement et al. 
1995). We will use this value for the mean absorp- 
tion of the LHB. It corresponds to approximately 
3.6% absorption of the \ keV emission. Other 
investigators have used greater columns; Sidher 
et al. (1996) used 6.6 x 10 18 cm' 2 (~ 5.9% ab- 
sorption), Kerp (1994) and Pietz et al. (1998) used 
1 x 10 19 cm -2 (~8.7% absorption), and Parmar 
et al. (1999) used 7 x 10 19 cm -2 (~ 43% absorp- 
tion). The last of these values is completely in- 
compatible with optical absorption line studies. 

Temperature: The model temperature of the 
LHB can be determined from the R2/R1 (or 
R2/R1L) ratio in a direction in which absorp- 
tion by the galactic ISM blocks \ keV emission 
from more distant sources. Figure 16 demon- 
strates the difficulty of this problem. Numer- 
ous strong features visible in the | keV band, 
mostly supernovae and super-bubbles, penetrate 
the galactic absorption to contaminate the local 
emission, and increase the R2/R1 ratio. One can 
calculate the hardness ratio in a large window to- 
wards the anticenter in between known sources 
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Fig. 16. — The ~ keY hardness ratio as a function of galactic longitude for the region within 5° of the 
galactic plane. Brighter, hotter, absorbed features that “leak” into the R12 band are labeled. “Scan” refers 
to region with poor exposure due to a missing scan. The horizontal bar overlaying the data shows the mean 
hardness ratio over the given ranges of galactic longitude. 


of contamination, and one finds that the temper- 
ature is ~ logT = 6.10 (R2/R1 — 1.057). The 
temperature of the LHB towards the galactic cen- 
ter does seem to be higher (Snowden et al. 1990b, 
2000), but the contamination by the galactic bulge 
“leaking” through the absorption does not allow 
a determination with this data set. 

One might think to avoid most sources of con- 
tamination by sampling just off of the galactic 
plane. The only region that is near the galac- 
tic plane and has a substantial absorbing column 
(10 21 < Nh < 2.0 x 10 21 cm -2 ) is the region of 
the Taurus-Auriga molecular clouds. Absorption 
line data place these clouds at — - 120 pc, so there 
is little room for extra emission components be- 
tween the LHB (radius ~ 100 pc in this direc- 
tion) and the absorbing clouds. A sample of the 
RASS of ~ 360 square degrees yields a R2/R1 = 
1.02 ± 0.01, or an equivalent Raymond k Smith 
model temperature of logT = 6.09 ± 0.01 

Thus, the current analysis produces an LHB 
temperature through broad-band ratios that is 
consistent with independent measures also made 
with broad band ratios. The significantly lower 
temperatures found by the spectral fitting of past 
analyses are due to the need for a much higher 
temperature TAE, when one assumes that the 
TAE is due to a single thermal component. 

5. Discussion 

Historical Background: The discovery of the 


soft X-ray emission that we have come to associate 
with the LHB occurred as the result of attempts to 
measure the soft extrapolation of the extragalactic 
power law (Bowyer et al. 1968; Bunner et al. 1969: 
Davidsen et al. 1972). The ubiquity of the £ keV 
emission was a mystery that was not to be solved 
by the discovery of instrumental effects or incor- 
rect modeling of the particle background, though 
these problems remain. (See, for example, the dis- 
cussion of the Long Term Enhancements and the 
earth as a “flaming comet” in Cox 1998.) The 
\ keV emission anti-correlation with the total neu- 
tral hydrogen column led to two main schools of 
thought. The first, the absorption model (Bowyer 
et al. 1968; Bunner et al. 1969; Davidsen et al. 
1972), was that the ubiquitous emission was extra- 
galactic and absorbed, though it was found that 
the absorption cross-sections could not produce 
the observed dynamic ranges (McCammon et al. 
1976; Long et al. 1976; Burrows et al. 1984) or 
band ratios (Bloch et al. 1986; Burrows 1989; Juda 
et al. 1991). The second, the displacement model 
(Sanders et al. 1977; Cox k Snowden 1986; Cox 
k Reynolds 1987; Snowden et al. 1990a), was that 
the emission was proportional to the unabsorbed 
path length, which was roughly anti-correlated to 
the total column. There was a third model, the 
embedded cloud model, Jakobsen k Kahn (1986), 
based somewhat on the McKee k Ostricker (1977) 
model of the ISM, that suffered from many of the 
same problems as did the absorption model (Bur- 
rows 1989). The advent of X-ray shadows (Snow- 
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den et al. 1991; Burrows k Mendenhall 1991), ob- 
served early in the ROSAT era, showed that both 
models were, to some extent, correct. 

It was clear that the soft X-ray background was 
composed of at least three components, the LHB, 
the EPL, and some component beyond the galac- 
tic absorption, which we have called the TAE. Al- 
though the Draco shadows (( = 89° - 91 °, b = 
38° - 40°) were the first glimpse of a hot galactic 
halo, it was not the long-theorized galactic corona 
(Spitzer 1956); the \ keV TAE was too patchy 
(Burrows k Mendenhall 1991; Snowden et al. 
1991, 1994a). It was suggested that the £ keV 
and | keV TAE were due to different components 
(Burrows k Mendenhall 1991) as the | keV emis- 
sion was much smoother than the | keV emission. 
Like the infamous Leiden-Groningen H I 21 cm 
data split that produced the artificial separation 
between high and intermediate velocity clouds, the 
RASS was split between the | keV and | keV 
bands and analysed by different groups. 

Summary of the Argument: The current work 
attempts to effect the obvious synthesis of the two 
bands. Our results have shown that the TAE is 
composed of at least two thermal components, 
perhaps more. The result is based on five argu- 
ments: 1) That the TAE of the R12 and R45 bands 
are not correlated, 2) That no single component 
can produce both the observed anti-correlation of 
emission and absorbing column and a consistent 
SED, 3) That the two-component models fit the 
SED better than any single-component model, 4) 
That single component models are not acceptable 
because they do not allow the LHB to have a tem- 
perature that is consistent with that measured by 
the R2/R1 ratio in the galactic plane where there 
is no “leaking” background emission, and 5) That 
the two component model is more consistent with 
detailed spectral fitting than are the single compo- 
nent models. The argument that the | keV emis- 
sion is much smoother than the ~ keV emission is 

4 

no more than suggestive. With the exception of 
such strong features as the Draco enhancement, 
a single-temperature model, in order to get the 
proper R45/R12, will necessarily scale down the 
structures seen in the R12 band. Since much of 
the emission in the R45 band is due to the EPL, 
the small variations due to the TAE would become 
difficult to detect in the R45 band. 

Although we posit two (or more) components 


to the TAE, the components are not neatly seg- 
regated into the \ keV and | keV bands; —22% 
of the | keV band flux (before absorption) is pro- 
duced by the soft component, while —12% of the 
£ keV band flux (also before absorption) is pro- 
duced by the hard component. Because the soft 
component does not have a significant trend with 
galactic longitude, global analyses that utilize only 
the R4-R7 bands, such as Wang (1997), will not 
be significantly affected by our analysis. 

Global Properties of the TAE: The mean surface 
brightness of the soft TAE component is — 816 
counts s -1 arcmin -2 or — 3.3 x 10 -8 ergs cm' 2 s“* 
sr -1 over the 0. 1-2.0 keV range; for the hard com- 
ponent, — 155 counts s -1 arcmin -2 or — 1.2 x 10“ 8 
ergs cm -2 s” 1 sr _1 over the 0. 1-2.0 keV range. 
Even away from Loop I (t — 330° ± 60°) and the 
Draco enhancement (( — 90°), the soft component 
is quite mottled, though the azimuthally averaged 
distribution is quite constant. The hard compo- 
nent, away from Loop I, is smooth (see Figure 9). 

These characteristics suggest that the hard 
component might be identified with the galactic 
corona or a galactic halo in hydrostatic equilib- 
rium (Wang 1997; Pietz et al. 1998), though that 
identification is by no means certain. If the hard 
component is the corona, with what is the soft 
component to be identified? 

The missing piece in halo studies is often the 
distance to the emission, and this study is no ex- 
ception. Several studies have set the £ keV ter- 
minus a quo. The ^ keV shadowing clouds in 
Draco have distances 2 > 180 pc (Lilienthal et al. 
1991). Benjamin et al. (1996) found distances of 
2 — 285 ± 75 pc to an infrared cirrus cloud in 
the Lockman hole that shadows the distant X-ray 
emission. The determination of the terminus ad 
quern is more difficult. The attempts to use high 
and intermediate velocity clouds are difficult as 
the disk absorption is much stronger, and there 
have been suggestions that the high velocity clouds 
may have X-ray emitting bow-shocks. A search for 
\ keV shadows due to high velocity complex M 
(1.5 < 2 < 4.4 kpc) was equivocal (Herbstmeier 
et al. 1995). A trivial upper limit is set by the 
shadowing of the extragalactic background by sev- 
eral galaxies (Barber & Warwick 1994; Snowden k 
Pietsch 1995; Barber et al. 1996; Cui et al. 1996). 
Distance determinations for the f keV emission 
by shadowing are even more difficult as high and 
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intermediate velocity clouds do not have sufficient 
column density to produce a measurable shadow. 

External Galaxies: Observations of external 

galaxies provide further clews to the distribution 
of the TAE. The - keY flux due to the TAE 
(both components) is ~ 750 counts s -1 arcmin -2 , 
while the | keV flux is ^ 40 counts s* 1 arcmin -2 . 
Thus, despite the absorption by our own galaxy, 
the closer external galaxies should have detectable 
halos. 

The X-ray halos of face-on galaxies are contam- 
inated by the point sources in the disk, but this 
contamination has diminished as resolution has 
improved. The best study to date is that of the 
closest face-on Milky Way analogue, M101, which 
was observed by the ROSAT PSPC (Snowden 
Pietsch 1995). If we use the relative D25 values for 
the two galaxies to scale to the equivalent of the 
galactocentric radius of the sun, a radius of ~ 8.'5, 
and if we assume that the entire M101 flux, after 
point source removal, is due to the halo, we find 
1570 counts s -1 arcmin -2 at - keV, while the flux 
at | keV at that radius is indistinguishable from 
the background. Significant \ keV flux (~ 350 
counts s -1 arcmin -2 ) extends to ~ 20'. (The scale 
length of the £ keV emission is ~ O.4.D25/2, and 
the scale length of the | keV emission would ap- 
pear to be similar.) Thus, both the strength and 
extent of the soft component of the TAE is con- 
sistent with that observed in Ml 01. The temper- 
ature determined from the R2/R1 ratio for the 
M101 halo is, however, cooler, logT = 5.86. 

Edge-on spiral galaxies for which extra-planar 
X-ray emission has been detected are often strong 
star-burst galaxies. Three non-starburst edge-on 
spiral galaxies have ROSAT - detected X-ray ha- 
los: NGC 891, NGC 4631, and NGC 4565. Of 
these, NGC 4565 has the lowest stellar formation 
rate (SFR, as measured by Lfir/D 2 5 ) and has 
no detected extended diffuse ionized gas (DIG). 
Vogler et al. (1996) detected some extraplanar dif- 
fuse X-ray emission, but were unable to determine 
a scale-height, and could only say that the tem- 
perature was consistent with a few times 10 6 K 
gas. NGC 4631, which has a higher SFR, is a dis- 
turbed galaxy that may be interacting (Weliachew 
et al. 1978; Rand et al. 1992; Rand & Stone 1996; 
Donahue et al. 1995). It has a strong nuclear star 
formation, and nuclear Ha filaments, which may 
be related to a £ keV emission region. The | keY 


emission has a HWHM z-extent of almost 3 kpc, 
a radial extent of about 3 kpc (0.25 £>25/2), and a 
temperature log Trs = 6.48 (Wang et al. 1995) or 
log Tmkl — 6.31 (Dahlem et al. 1998). NGC 891, 
which has an even higher SFR, has a more uni- 
formly extended X-ray halo; the | keV emission 
has a HWHM 2-extent of 2.1 kpc, a radial scale 
length of 4.1 kpc (0.3 D25/2), and a temperature 
of log T = 6.54 (Bregman &: Houk 1997). 

The only extragalactic detection (in an edge-on 
spiral) of a component that might be similar to the 
soft component of the TAE is that by Wang et al. 
(1995) in NGC 4631, where they detect a com- 
ponent with logT < 5.78. Dahlem et al. (1998) 
also detect a component with logT = 5.76, but 
Read et al. (1997) do not. Such a soft compo- 
nent could not be detected in NGC 891 due to 
the higher Galactic absorption. Lack of detection 
of the soft TAE component in external galaxies is 
not surprising; it is subject to strong absorption 
by both the Milky Way and the host galaxy. If the 
soft component has a lower scale-height than the 
hard component, then not only will the absorption 
by the host galaxy be important, but the ROSAT 
resolution would not allow separation of the soft 
TAE from disk emission. 

Pure Speculation: Both NGC 4631 and NGC 
891 have SFRs, based upon Lfir/D f 5 , that are 
similar to that of the Milky Way (Rand 1998). 
Both NGC 4631 and NGC 891 have halo temper- 
atures that are consistent with the logT = 6.42 
that we have found for the hot component of the 
TAE. Given the similarities of halo temperature, 
and that the extent of the extra-planar diffuse ion- 
ized gas is roughly dependent upon the SFR, it 
seems likely that the smooth hard component of 
the TAE has z-scale heights ~ 2 kpc, somewhat 
smaller than the 4 kpc exponential scale height 
fitted by Pietz et al. (1998). Both M101 and 
NGC 891 imply that the hard component should 
have a radial scale length ~ O.35D05/2. or ~ 4 
kpc, which is much smaller than the 15 kpc scale 
length fitted by Pietz et al. (1998), but implies a 
much more disk-like distribution than the model 
of Wang (1997). 

The | keV emission in the Milky Way is quite 
symmetric, the mean flux in the south being very 
similar to that of the north, the emission is quite 
symmetric in NGC 891, and the asymmetry of 
NGC 4631 is thought to be an inclination effect 
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(Wang et al. 1995). With a disk-like distribution, 
this symmetry would argue that the hard compo- 
nent is hot gas in hydrostatic equilibrium, where 
the symmetry across the disk is enforced by the 
gravitational potential. 

The soft component, however, is very differ- 
ent. The soft component is patchy, and there is 
a distinct north-south asymmetry in the | keV 
TAE; the northern hemisphere is about 1.9 times 
stronger than the southern hemisphere (Snowden 
et al. 1998). and the temperature of the soft com- 
ponent in the southern hemisphere may be lower 
than that of the north (Snowden et al. 1998) as 
is implied by Figure 11 of Snowden et al. (2000). 
These properties suggest that the soft component 
is the result of local phenomena. Perhaps the soft 
component is composed of the plumes of hot gas 
ejected from galactic chimneys. Snowden et al. 
(1998) suggested that the Draco enhancement was 
the extraplanar breakout of the Cygnus superbub- 
ble, though the intervening mid-latitude absorp- 
tion makes testing this model difficult. (It is in- 
teresting to note that such a model suggests that 
the Draco emission occurs to 2 ~ 2 kpc.) How- 
ever, as the j keV TAE is not restricted to lower 
latitudes any local phenomenon must be fairly uni- 
formly distributed or must send gas to fairly high 
z. Shelton (1998) has suggested that supernovae 
remnants in the halo might dominate the soft X- 
ray emission and produce the patchy structure. 

As has been pointed out by Gehrels & Williams 
(1993) in the context of stellar coronae, the cooling 
curve for a hot thin plasma has two local minima 
at temperatures that would emit in the soft X- 
ray range. Of course, the exact temperatures of 
those minima (logT = 5.9 and 6.6) do not agree 
particularly well with the observed temperatures 
(logT = 6.06 and 6.42). Reference to Schmutzler 
& Tscharnuter (1993) shows that the location of 
the minima of the cooling curve are sensitive to the 
metallicity and previous ionization history of the 
gas. It is only suggestive that the minima of the 
cooling curve may correspond to the temperatures 
of the dominant components. In this case, the 
ubiquity of the soft component suggests that the 
soft component may represent some part of the 
hard component at a more advanced age. 

Future Work: Determining the true thermal 
structure of the galactic halo will be difficult. 
Next generation X-ray instruments will have lower 


background contamination, but typically will have 
small fields-of-view, so using angular variations 
caused by absorption to separate local and distant 
components will be more difficult. Although spec- 
tral resolution will be good at f keV, the spectral 
resolution below the carbon edge, where the LHB 
and the soft component of the TAE has most of 
their emission, remains poor. Experience with the 
DXS has shown that given high spectral resolu- 
tion, understanding the spectrum of the LHB is 
still difficult, and the parameter space allowed by 
non-equilibrium plasma codes is large. 

Attention will also need to be given to incorpo- 
rating the X-ray TAE into our understanding of 
the halo optical/UV emission/absorption line ra- 
tios. Many line ratios are adequately explained by 
photoionization due to massive disk stars, whose 
ionizing radiation is absorbed and hardened by 
passage through the gaseous disk (see, for exam- 
ple, the review by Rand 1998). However, there are 
still outstanding difficulties, such as the run of line 
ratios with z-height, in which the X-ray halo may 
play a part. 

We would like to thank Robin Shelton, Ian 
George, Richard Mushotzky, Peter Teuben, Richard 
J. Rand, and Sylvain Veilleux for several interest- 
ing discussions. This study was funded in part by 
the NASA ADP program. 
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